Upper and lower critical fields in NdFeAs(0,F) single crystals : a study by Hall probe 

magnetization and specific heat 
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The upper and lower critical fields have been deduced from specific heat and Hall probe magne- 
tization measurements in non-optimally doped NdFeAs(0,F) single crystals (T c ~ 32 — 35K). The 
anisoptropy of the penetration depth (I\) is temperature independent and on the order of 4.0± 1.5. 
Similarly specific heat data lead an anisotropy of the coherence lenght ~ 5.5 ± 1.5 close to T c . 
Our results suggest the presence of rather large thermal fluctuations and to the existence of a vortex 
liquid phase over a broad temperature range (~ 5K large at 2T). 

PACS numbers: 74.60.Ec, 74.60.Ge 



The recent discovery of superconductivity at unusually 
high temperature (up to 56K) in rare earth iron oxypnic- 
tides 0, |2] has been the focus of a tremendous number 
of theoretical and experimental works in the past few 
months. The possible coexistence of superconductivity 
with a complex magnetic structure [|[ makes this system 
particularly fascinating and a non conventionnal pairing 
mechanism (associated with multi-gap superconductivity 
in electron and hole pockets) has been suggested by dif- 
ferent groups In this context, it is very important to 
have a precise determination of the temperature depen- 
dence of the critical fields and corresponding anisotropics. 

We will focus on the superconducting properties of the 
Nd(0,F)FeAs compound [5]. The aim of this paper is to 
present combined Hall probe magnetization and specific 
heat (C p ) measurements performed on Nd doped single 
crystals. A superconducting anomaly was clearly visi- 
ble in C p around 35 K and, as previously observed in 
high T c oxides, this anomaly is symmetric, broad and 
rapidly collapses as the magnetic field is increased. On 
the other hand, the lower critical field presents a clas- 
sical temperature dependence for both ff ||c and H a \\ab 
(with H^(0) ~ 120 ± 30G and ~ 40 ± 10 G) and 
the anisotropy of the penetration depth remains on the 
order of 4 on the entire temperature range. A similar 
anisotropy is obtained for the upper critical field from 
the Cp measurements. 

Nd(0,F)FeAs samples have been synthetized at high 
pressure in a cubic, multianvil apparatus. More de- 
tails of the synthesis are given elsewhere [6| . Both Hall 
probe magnetization and C p measurements have been 
performed on the same platclct-likc single crystals ex- 
tracted from the polycrystalline batch [with dimensions 
: - 100 x 100 x 30 ^m 3 (sample 1), - 120 x 80 x 30 
(sample 2) and - 170 x 100 x 50 (sample 3)]. 

The local magnetization of the platelet has been mea- 
sured by placing the sample on the top of an array of 



miniature Hall probes of dimensions 4 x 4 or 8 x 8/im 2 . 
The AC transmittivity (T' ac , related to the local suscep- 
tibility) has been measured by applying a small (~ 1G) 
AC field and recording the corresponding response of the 
probe as a function of the temperature (see Fig. la). T' ac is 
obtained by subtracting the response in the normal state 
from the data and is rescaled to — 1 at T = 4.2K (due to 
the small but non zero distance between the probe and 
the sample surface about 10% of the external field is still 
picked up by the probe in the superconducting state). 
As shown in Fig. la, the transmittivity measured in the 
center of the sample first presents a small paramagnetic 
increase starting at ~ 37K followed by a sharp diamag- 
netic jump at T ~ 34K (sample 1). This paramagnetic 
bump (also observed in zero field cooled DC measure- 
ments) reflects a non homogeneous distribution of the 
field between 34 and 37 K 0] . However the presence of a 
clear anomaly in the specific heat emphasizes the overall 
good quality of the platelets. A similar bump has been 
observed in sample 2 between 34. 5K and 37. 5K and in 
sample 3 between 36. 5K and 37.5 K. Note also that pre- 
vious global magnetization measurements @ yield to an 
onset of the diamagnetic response of the polycrystalline 
batch around T ~ 51K in agreement with a resistivity T c 
on the order of 50K for optimally doped samples Q . We 
did not find any sign of superconductivity around 51K 
in our platelets indicating that the diamagnetic response 
observed above 40K in global measurements is due to a 
distribution of T c values (i.e. F content) within the poly- 
cristalline batch. 

Specific heat measurements have been performed on 
the same crystallites (sample 1 and 2) using an AC tech- 
nique. This high sensitivity technique (typically 1 part 
in 10 4 ) is very well adapted to measure C p of very small 
samples. Heat was supplied to the sample at a frequency 
lu = 10 Hz by a light emitting diode via an optical fiber. 
The temperature oscillation have been recorded with a 



2 




25 30 35 40 



T(K) 

FIG. 1: (color online) (a) Temperature dependence of the AC 
transmittivity and specific heat for the indicated field val- 
ues (i? a ||c). (b) Temperature dependence of the specific heat 
in a Nd(F,0)FeAs platelet (sample 1) for i/ a ||c (closed black 
symbols) and H a \\ab (open red symbols) for the indicated field 
values. The inset shows sample 1 mounted on the thermocou- 
ple cross, the width of the legs is ~ 50/xm as well as sample 3 
(bottom left corner). 



thermocouple which has been calibrated from measure- 
ments on ultra pure silicon. The superconducting con- 
tribution to the specific heat (AC P ) has been obtained 
by subtracting the curve at fioH a = 7T (for H a \\c) from 
the curves obtained for lower fields as well as a (H a /T) 2 
contribution to account for the presence of a magnetic 
background . As shown in Fig. lb (sample 1), there is 
no sharp discontinuity at T c and the anomaly has the 
shape of a rounded peak with a height on the order of 
a few 10~ 3 of the total specific heat. The other strik- 
ing behaviour is the influence of the magnetic field. As 
previously obsersed in cuprates [11] , the smearing of the 
anomaly with field is unexpectedly strong with an al- 
most field independent onset. The maximum of C p is 
pushed downwards in temperature by the field with a 
characteristic average slope on the order of —2 T/K for 
H a \\c and —7 T/K for H a \\ab) [compared to resp. —1.5 



T/K and -9 T/K in optimally doped YBaCuO, note 
that the shift is possibly non linear for small magnetic 
field, see solid line in Fig. 2 and discussion below]. As in 
cuprates, the shape of the anomaly precludes any pre- 
cise determination of H C 2- However, assuming that H C 2 
is close to the maximun of C p {H max ) and taking the 
corresponding linear temperature dependence, we hence 
get H c2 (0) ~ 0.7 x dH max /dT x T c ~ 40T and ~ 170 T 
for -ff a ||c and H a \\ab respectivevely. The corresponding 
anisotropy of the coherence length = £ c /£ a b = T Hc2 
is hence on the order of H^(0)/H^ 2 (0) ~ 4. The uncer- 
tainty on the H C 2 values is rather large but can also be 
estimated from the total shape of the anomaly writting 
: AC(r^ x H a \\ab) ~ AC(H a \\c) consistently leading to 
~ 5.5 ± 1.5. This value is also in good agreement with 
the one previously reported by Welp et al. [13[ and with 
recent band structure calculations [14j ]. 

The influence of H a on the C p anomaly and its overall 
shape is strikingly different from the one observed is clas- 
sical superconductors or even in the (K,Ba)BiC"3 system 
(T c ~ Z2K) @ or MgB 2 (T c ~ 39 K) [lo| of similar T c val- 
ues. It is tempting to attribute these deviations to fluc- 
tuation effects. To reinforce this idea, it is worth noting 
that the onset of the diamanetic response well coincides 
with the inflection point of the C p anomaly for H a = 
but is pushed towards substantially lower temperature in 
field (~ 5K lower than the C p shift for t i H a = 2T). As 
T' ac is sensitive to pinning, this result suggests that the ir- 
reversibility line is well separated from the superconduct- 
ing transition for H a ^ (see inset of Fig. 2 for for H a \\c, 
a similar result - not shown - is obtained for _ff a ||a6) sup- 
porting the idea of an extended vortex liquid phase. The 
importance of thermal fluctuations can be quantified by 
the Ginzburg number [l2| Gi = (fcBT , c /e £ c ) 2 /8 where e 
(= (<&o/47rA a b) 2 ) is the line tension of the vortex matter. 
Large A values (see below) combined with small £ values 
hence lead here to eo£c ~ 200 K (i.e. similar to cuprates) 
compared to ~ 1000 K in (K,Ba)BiC>3 and even ~ 10 4 in 
MgB2. Although smaller than in cuprates due to smaller 
T c values, thermal fluctuations are hence expected to play 
a role in this system possibly leading to the melting of the 
vortex solid (Gi ~ 3 x 10~ 3 — 10 -2 ). A similar conclusion 
has been drawn by Welp et al. [13j |. Note that, as ex- 
pected for vortex melting, the irreversibility line clearly 
presents a positive curvature (Hi rr oc (1 — T/T c ) a with 
a ~ 1.5, for a review see [12j). Similarly, thermal fluctu- 
ations may also induce an upward curvature in the H C 2 
line [IH but eventhough our data suggest the existence of 
such a curvature (see also [Hf), we cannot unambiguously 
conclude given the uncertanties and limited temperature 
range of the C p data. 

The first penetration field H p has been deduced by 
measuring the remanent field (B rem ) in the sample af- 
ter applying an external field H a and sweeping the field 
back to zero. For H a < H p no vortices penetrate the 
sample and the remanent field remains equal to zero. 
H a is progressively increased until a finite remanent field 
is obtained as vortices remain pinned in the sample for 
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FIG. 2: (a) Temperature dependence of the upper critical field 
(maximum of the specific heat anomaly, see Fig. fa) for // a ||c 
(circles) and H a \\ab (squares). In the inset : Remanent field 
(B rern ) as a function of the applied field (J,oH a (see text for 
details) for the indicated temperatures and field orientation. 
The solid lines are (H a — H p ) 2 fits to the data, (b) : compar- 
ison between H C 2 (deduced from C p , closed symbols) and the 
onset of the diamagnetic response Hi rr (open symbols) for 
H a \\c suggesting the existence of a broad vortex liquid phase. 



H a > Hp (B rem increasing as {H a — H p ) 2 , see solid line 
in the inset of Fig. la). Theoretically speaking, H p is 
expected to depend on the position of the probe, being 
larger in the center of the sample as vortices first remain 
pinned close to sample edges. However, experimentally, 
a non-zero B rem value could be detected almost simul- 
taneously on all the probes due to the non zero distance 
between the sample and the probe. Note that we did not 
observe any significant change in the B rern (H a ) curve up 
to H a ~ 3H p as the magnetic field was tilted away from 
the c-axis (with an angle 6jj a < 60 — 70°). Even though, 
the angle between the internal field {9h) and the c-axis 




T/T 



FIG. 3: Temperature dependence of the lower critical field 
Hd along the c direction (solid squares for sample 3, solid 
lozanges for sample 1, rescaled by a factor of 1.5) and ab plane 
for sample 3 (solid circles). The solid lines are the standard 
behaviours expected from the the BCS theory. The corre- 
sponding temperature dependence of the anisotropy (r^ = 
1.3 x r_ff c i) is also reported (open squares). The error bars 
correspond to the uncertainties on the determination of H p , 
the whole curve might be shifted ±1.5 due to the uncertainty 
of the demagnetization factor. In all cases, Tx remains flat in 
T. The Fj values deduced from C p (present work : open circle 
and ref [13] : open lozange) and transport : open triangles 
(from ref. [26]) measurements are also displayed. 



is reduced by the demagnetization effects, this surprizing 
behaviour suggests that the induction {B) in the sample 
remains perpendicular to the platelets up to large Or val- 
ues (see inset of Fig. la for sample 1 at 9n a = 7° and 55°). 
B rem finally decreases with 8 for H a > 2>H p (and/or for 
9n a > 70°), indicating that the vortices finally aligns on 
the applied field for large H a (and/or 6># a ) values. 

In samples with rectangular cross sections, flux lines 
partially penetrate into the sample through the sharp 
corners even for H a < H p pjl, [l7| but remain " pinned" 
at the sample equator. The magnetization at H a = H p 
is then larger than H c i and the standard " elliptical" cor- 
rection for H c i (= H p /(1 — N) where N is the demag- 
netization factor) can not be used anymore. Following 
[l6|, in presence of geometrical barriers, H p is related to 
H c i through H c \ w H p /tanh(^ad/2w) where a varies 
from 0.36 in strips to 0.67 in disks (2w and d being the 
sample width and thickness, respectively) .Taking an av- 
erage a value ~ 0.5 [HI we hence got [19| : H c cl {Q) ~ 
110 ± 20G in sample 1, ~ 100 ± 20G in sample 2 and 
H^(0) - 140 ± 30G in sample 3. To obtain H p in the 
ab- plane, sample 3 has been rotated by 90° putting the 
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probe on the small side of sample (i.e. applying H a along 
the long direction of the sample so that H^ b {Q) ~ H^). 

We hence obtained Hgf ~ 40 ± 10G (in sample 3). 

The temperature dependence of H c \ for both H a \\c 
(sample sample 1 and 3; the data for sample 1 have been 
rescaled by a factor 1.5) and H a \\ab (sample 3) is dis- 
played in Fig. 3. As shown, clearly flattens off at low 
temperatures (down to 1.6K for sample 1) indicating that 
the gap is fully open in our Nd(0,Fe)AsFe single crystals. 
This result is in striking contrast with H p measurements 
by [1H which suggested that H c \ varies linearly at low 
temperature in La(0,F)FeAs hence suggesting the exis- 
tence of gapless superconductivity. On the other hand, 
a very similar temperature dependence (and iJ c i(0) val- 
ues) has been recently obtained by Okazaki et al. in Pr 
doped samples [2^ |. 

The lower critical field is related to the penetration 
depth through : H c cX = $§/ ' (ATr\ 2 ab )(Ln(K) + c(k)) where 
k = Xab/£,ab (£ being the coherence length) and c(k) a 
k dependent function tending towards ^0.5 for large k 
values. Taking iJ c c 2 (0) - 40T, one gets A q6 (0) - 270 ± 
40 nm [2(J. Those A values are consistent with those 
previous obtained in oxypnictides. Indeed, slightly lower 
Xab values (~ 200 nm) have been reported in both Sm and 
Nd doped samples with higher T c s (~ 50K, from either 
torque |26| . muons relaxation [24[ or plasma frequency 
[25[ | measurements) whereas A ~ 360 — 390nm have been 
obtained in La doped samples with lower T c values (~ 
23K, 

For H\\ab, H$ = $§/(47rA a6 A c )(Ln(/c*) + c(k*)) with 
k* = A c /U, and H$ = $2/(2< c £ a6 ) ~ 170 T, leading 
to A c ~ 1200 ± 200 nm and a corresponding anisotropy 



(r^) on the order of 4.0 (±1.5), almost temperature in- 
dependent down to 1.6K (see Fig. 3, sample 3, T\ — 
T Hcl x (1 + Ln(T x )/(Ln( K ) + 0.5)) ~ 1.3 X T Hcl ). Note 
that our data are in striking contrast with those de- 
duced from torque measurements in Sm doped samples 
f26| which led to strongly increasing T\ values, exceed- 
ing 30 at low temperature. On the other hand, our T\ 
consistently tends towards for T — > T c deduced ei- 
ther from C p (present work and [l3j]) or transport data 
[27l ]. The difference between the anisotropics measured 
at low (r\) and high (r^) fields remain within our error 
bars and we hence cannot comment on a possible similar- 
ity with MgB2 in which the coexistence of 2 gaps leads 
to a strongly field and temperature dependence of the 
anisotropy [28j |. 

In conclusion, specific heat and Hall probe magnetiza- 
tion data suggest than r ? ~ T x ~ 5 ± 2 in Nd(F,0)FeAs 
single crystals. Whereas the irreversibility is rapidly 
shifting towards lower temperature for increasing H a , the 
C p anomaly collapses and shows only a minor shift sug- 
gesting the presence of strong thermal fluctuations pos- 
sibly leading to the melting of the vortex solid. 
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